import numpy as np
import matplotlib.pyplot as plt

mu = 1  # 均值
sigma = 3  # 标准差
num_samples = 1000000  # 随机数的数量

rand_data = np.random.gamma(shape=mu, scale=sigma, size=num_samples)  # 生成伽马分布随机数
plt.hist(rand_data, bins=100)  # 绘制直方图
plt.show()  # 显示图形